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Abstract — In this paper we show how to correct 
and improve the Belief Interacting Multiple Model fil- 
ter (BIMM) proposed in 2009 by Nassreddine et al. 
for tracking maneuvering targets. Our improved al- 
gorithm, called PCR-BIMM is based on results devel- 
oped in DSmT (Dezert-Smarandache Theory) frame- 
work and concerns two main steps of BIMM: 1 ) the 
update of the basic belief assignment of modes which is 
done by the Proportional Conflict Redistribution Ride 
no. 5 rather than Smets’ rule (conjunctive rule); 2) 
the global target state estimation which is obtained from 
the DSmP probabilistic transformation rather than the 
commonly used Pignistic transformation. Monte- Carlo 
simulation results are presented to show the perfor- 
mances of this PCR-BIMM filter with respect to clas- 
sical IMM and BIMM filters obtained on a very simple 
maneuvering target tracking scenario. 

Keywords: Tracking, Maneuvering target, IMM, 
BIMM, DSmT. 

1 Introduction 

In Fusion 2009 international conference, Nassreddine, 
Abdallah, and Denceux [13] have proposed an inter- 
esting idea to extend the classical Interacting Multiple 
Models (IMM) filter with belief function theory in order 
to deal with an unknown and variant motion models. 
Their algorithm is based on the classical/historical be- 
lief function theory developed by Shafer in 1976 [14], 
known as Dempster-Shafer Theory (DST) and requires 
both Smets’ rule, i.e. the conjunctive fusion rule equiv- 
alent to the non normalized Dempster’s rule, and the 
probabilistic pignistic transformation. This algorithm 
is called Belief Interacting Multiple Model algorithm 
(BIMM). According to authors results, BIMM algo- 
rithm outperforms classical IMM algorithm at least in 
the vehicle localization problem studied in their works. 
These appealing results and the possible extension of 
IMM in belief function theory framework motivates our 
interest to analyze and evaluate this new BIMM filter. 



A deep analysis of the paper yields to the following 
comments: 

1. The derivation of the predicted prior basic belief 
assignment of modes in Step 1 of BIMM algorithm 
was clearly wrong in [13] as proved in the sequel. 
This mistake implies a serious doubt on the validity 
of the results presented in [13]. 

2. The simulations results presented in [13] cannot 
be verified precisely, nor reproduced, because some 
settings parameters (like a, discounting factors) re- 
quired for the BIMM filter have not be provided 
by the authors and the essential step 9 of the algo- 
rithm was not detailed enough. 

3. It is known (see Chapter 1 of [15] Vol. 3) that 
the conjunctive rule does not perform efficiently in 
a sequential fusion process because the empty set 
is an absorbing element for the conjunctive fusion 
rule. Therefore, in order to implement successfully 
the BIMM filter, some ad-hoc numerical techniques 
are necessary (or some extra normalization steps) 
in the BIMM algorithm in order to prevent the 
mass of belief committed to empty set to become 
close to one and make Smets’ rule responding to 
new information. This serious problem has unfor- 
tunately not been discussed in [13]. 

From the theoretical point of view, it is quite surpris- 
ing that one gets better performances with the BIMM 
(which proceeds with less specific information since it 
deals with non Bayesian basic belief assignments) than 
with the classical Bayesian IMM filter (which deals with 
more specific information, i.e. with Bayesian basic be- 
lief assignments). The first purpose of this work is 
to verify if the conclusions given in [13] are valid on 
a very simple reproducing maneuvering target track- 
ing scenario. We want also to see if a more justified 
Belief-based IMM algorithm can be developed to im- 
prove the BIMM algorithm and to evaluate it to get 
a fair comparison of its performance with respect to 
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classical IMM filter. The improvement of the BIMM 
algorithm we propose in this paper is based on ad- 
vanced theoretical results obtained in the development 
of Dezert-Smarandache Theory (DSmT) of information 
fusion [15]. This paper is organized as follows: After a 
brief recall of classical (fixed structure) IMM algorithm 
given in section 2, one presents in section 3 the Belief 
IMM algorithm and its flaws. Motivations for the im- 
provement of the BIMM filter is presented in section 
4 with the presentation of the main steps of our new 
algorithm called PCR-BIMM filter (Proportional Con- 
flict Redistribution-based BIMM). In section 5, we ex- 
amine the performances of the IMM, and PCR-BIMM 
on a very simple tracking scenario through Monte- Carlo 
results. Conclusions and perspectives for further inves- 
tigations are given in section 6. 

2 Classical IMM algorithm 

The IMM filter is one of the most used algorithm for 
tracking maneuvering targets and was developed origi- 
nally by Henk Blom in eighties [5, 6, 2], The IMM filter 
is a recursive filter with a low complexity and has been 
proved very efficient in many real-data tracking appli- 
cations [4] and many extensions of IMM have been de- 
veloped since its original publication for dealing with 
nrultitarget-multisensor case, cluttered environments, 
etc, see [12] for a good survey of Multiple Models tech- 
niques. The classical IMM algorithm considers a hybrid 
Multiple Models (MM) system which obeys one of a fi- 
nite number r of dynamic models Mi, i = 1 , ,r and 
estimates the posterior mode probabilities from their 
prior probabilities and target measurements (Bayesian 
framework) . Its specificity is that IMM mixes hypothe- 
ses with depth 1 only at the start of each cycle and thus 
has a low complexity of order 0(r), while providing 
same performances as the more effective Generalized 
Pseudo-Bayesian estimator of order 2. We briefly recall 
the principle of classical IMM filter, see [3, 4] for more 
details with examples. A hybrid MM system is charac- 
terized by two state variables: 1) the base-state variable 
x(fc) of dimension n x including the position, velocity, 
etc. of the target, and 2) a modal-state Mj(k) belong- 
ing to a known finite set M r {k) = {Mi{k),i = 1, . . . , r} 
of r possible dynamic models for the target during its 
motion. For simplicity of presentation, we consider only 
a fixed-structure IMM, i.e. A4 r (k) = A4 r is invariant 
with time. Variable-structure IMM is possible and has 
been introduced by Xiao-Rong Li in [10, 11]. The hy- 
brid system is described by the equations 1 

x(k ) = F [M(k)]x{k — 1) + v[fc — 1 ,M(k)\ 

z (k) = H [M (k)]x(k) + w[/c, M (fc)] 

where M(k) is the mode in effect during the sampling 
period ending at time k belonging in A4 r . x(k ) and 

1 For simplicity, we assume here linear systems. 



z(k) are the target state and observation vectors. The 
set of all available measurements up to k is denoted 
Z k . F[M(k)] and H [M(k)\ are known matrices de- 

pending on the dynamic model M(k). The statistics 
of the process and observation noises v[fc — 1 ,M(k)] 
and w[/c, M{k)\ can differ from mode to mode. Usually 
one considers v[fc — 1 ,M(k) = Mj] ~ J\T(vj,Qj) and 
w[fc, M(k) = Mj] Af(wj,Rj) with known covariance 
matrices Qj and Rj respectively. The Mode jump pro- 
cess is modeled as a Makov chain with known a priori 
probabilities P{M(0) = Mj} = jij(k = 0) and known 
transition probabilities P{M(k ) = Mj\M(k — 1) = 
Mi} = nij. A cycle of the classical IMM algorithm 
(k — 1) i-> k consists in the following steps: 

• Step 0 (Initialization at k = 0): Definition of dy- 
namic and observation matrices, choice of process and 
observation noise levels, sampling period, initialization 
of the filters adapted to each mode, choice of the prior 
mode probabilities Pj and the transition probability 
matrix P t = [ir.ij = P{Mj(k)\Mj(k — 1)}]' assumed 
known and time- invariant. 

• Step 1 (Interaction- mixing (j = 1, . . . , r)): Mixing 
of the previous cycle mode-conditioned state estimates 
Xj(fc — l|fc — 1) and covariance, using the mixing prob- 
abilities Hi\j(k — l|fc — 1), to initialize the current cycle 
of each mode-conditioned filter x°(fc — l|fe — 1). This is 
done by 

r 

Xj (k — 1 1 A: — 1) = y^ Mi |.(fc-l|fc-l)x,(fc-l|fc-l) (1) 

i=l 



P° j (k-l\k-l)=Y,n i]j (k-l\k-l){P i (k-l\k-l) 

i= 1 

- [x-i(k - 1| k - 1) — x°j(k — 1| k — 1)]- 

[xi(k - 1| k - 1) - x°j(k - 1| k - 1)]'} (2) 

where the elements /iju(fc— l|/c — 1) of the mixing prob- 
ability (vertical) vector p fc _ 1 | fc _ 1 (.| Mj(k)) = [fJ>i\j(k — 
l\k — l),i = 1 , . . . r]' are calculated by 

m\j{k - 1| k - 1) 4 P{Mi(k - 1)1 Mj(k), Z k ~ 1 } 

= KijlMjk - 1 ) 

/u ( k ) 

with 

r 

(k) ^ P{M } (k) | z fc_1 } = £ t ( k - 1) (4) 

i = 1 

The equation (4) can be written more concisely as: 

Mfc (•) = p * • Mfe-i(-) (5) 

where P' t = [nij] and p fc _ 1 (.) represents the (vertical) 
vector of prior probability of modes, i.e. 

Vk-i(-) = [P(Mi(k—l)\Z k ~ 1 ] 1 = [ni(k—l) ...Mr(fc-l)]' 
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and pu k (.) represents the (vertical) vector of predicted 
prior probability of inodes 

%(•) = [P(Mj(k) IZ*- 1 ]' = fa (k) ...fa(k)Y 

• Step 2 (Mode conditioned filter): From prior 

mixed statistics x 6 (fc — l|fc — 1) and P°(fc— l\k— 1) and 
the target measurement z (k), one calculates itj(k\k) 
and Pj(k\k) for each possible mode in effect (r filters 
running in parallel) by a specific filter matched to mode 
Mj , typically a Kalman filter if the dynamic and obser- 
vation system are linear, or Extended Kalman Filter 
(EKF) to deal with linear or non linear equations, or 
any other sophisticated filters if necessary for dealing 
for example with miss-detections and false alarms [3]. 
The likelihood A j(k) of the filter j is assumed to be 
Gaussian with 



A j(k) = 



(2n) n */ 2 ^/\SjW)\ 



exp 



(k)Sj 1 (k)ij(k) 



where z j(k) = z(k) — Zj{k\k — 1) is the innovation and 
S j(k) is the covariance of the innovation provided by 
the filter j. 

• Step 3 (Mode probability update): The probability 
/ij(/c) of each mode j for j = 1, . . . , r is calculated by 



H {k) = P{Mj(k)\Z k } = Aj(k)fa(k)/ ^^Ai(k)fa(k) 

i = 1 

(7) 

• Step 4 (Global estimation for output purpose): 
The global estimate x(fc|fc) and the covariance of esti- 
mation error P(fc|fc) are given by: 



state of classical IMM filter is replaced with the pignis- 
tic averaging of the mode-conditioned state estimates. 
This new extension of IMM filter was called BIMM 
(Belief-based IMM) since it uses belief function theory 
to represent the uncertainty in the switches between 
the modes. This section presents succinctly the prin- 
ciple of the BIMM filter. We justify also our motiva- 
tion for developing a new Belief-based IMM algorithm. 
The steps of BIMM are actually very close to the steps 
of classical IMM, except that predicted and updated 
mode probabilities are estimated from pignistic proba- 
bilities derived from a basic belief assignment updated 
with the conjunctive rule of combination. The main 
changes of BIMM concern the Step 1 and the Step 3 
of IMM algorithm. The frame of discernment chosen 
in BIMM coincides with the set of possible models, i.e. 
0(fc) = M r {k) = {Mi(k),i = l,...,r}. Instead of 
computing recursively the mixed Pi\j(.) and updated 
Pj(.) probabilities with eqs. (3) and (4) as done with 
the classical IMM, one deals with bba’s defined on the 
power-set 2 e of the frame of discernment. Mathemati- 
cally, a normal bba m (.) is defined 4 as a mapping from 
2 e i— > [0, 1] such that m(0) = 0 and ^T 4g2 © m (^) = 1- 
A is a focal element of m(.) if m(A) > 0. Any discrete 
probability measure can be interpreted as a special be- 
lief function, called Bayesian belief [14] whose focal el- 
ements are singletons of 2 e . Any belief function with 
a bba m(.) can be approximated into subjective prob- 
ability measure thanks to the pignistic transformation 
[17] defined for all Mj G 0(fc) by 



BetP{Mi } 



E 

A£2 e \AnMi=Mi 



1 

Ri 



m(A) 

1 — m(0) 



(10) 



it(k\k) = ^2n j (k)± j (k\k) (8) 

3 = i 
r 

P(fc|fc) = ^/x,(fc){P i (fc|fc) 

i= i 

- [*#’ \k) - x(fc|fc)] • [x 3 {k\ k) - x(fc|fc)]'} (9) 

3 Belief-based IMM algorithm 

In 2009, Nassreddine et al. have proposed in [13] 
an extension of classical IMM filter in the framework 
of Dempster-Shafer Theory (DST) [14] for dealing with 
an unknown and variant motion models. The idea was 
to select a set of candidate models 2 , and then esti- 
mate a current basic belief assignment (bba) defined 
on the power-set of this set of models based on the fu- 
sion of bba’s built from measurement likelihoods with 
the predicted bba of the models using Smets’ rule 3 de- 
noted ©. From the result of Smets’ fusion, the mixed 

2 Corresponding to the so-called frame of discernment and usu- 
ally denoted © in DST. 

3 Smet’s rule is nothing but the non normalized Dempster’s 
rule of combination, i.e. the conjunctive rule. 



where |A| is the cardinality of A. 

The steps of BIMM proposed in [13] are 5 : 

• Step 0 (Initialization at k = 0): Definition of dy- 
namic and observation matrices, choice of process and 
observation noise levels, sampling period, initialization 
of the filters adapted to each mode. The prior proba- 
bilities of modes { Pj = P{M{ 0) = Mj},j = 1 ,...,r} 
used in IMM, are replaced 6 by the vacuous belief as- 
signment m(0(fc = 0) = Mi U AI-i U . . . U M r ) = 1. The 
probability transition matrix P( = [7Tj j\ is replaced by 
a bba transition matrix' M t = [m,;j] having a very sim- 
ple structure defined by the r implication rules: ” Rp. if 
AI(k) = Afj(fc) then M(k + 1) = Mj(fc + 1)” with known 
belief coefficients /% G [0,1] for i = 1,2 ,...,r with 
/3j = m(Mj(fc + l)|Mj(fc)) and 1 — /3j = m(0(fc + 1) = 
Mi U . . . U M r {k + l)|Mj(fc)). 

4 We use boldface letters to denote vectors or matrices. 

5 We use a more classical notation generally adopted in the 
tracking community. 

6 Note that this initialization can also be done by taking 
m(Mj(k = 0) = Pj as well if one considers that prior proba- 
bilities of modes is accurate enough. 

7 Called switching mass function in [13]. 



323 




Advances and Applications of DSmT for Information Fusion. Collected Works. Volume 4 



• Step 1 (Interaction- mixing): The mixing probability 

— 1 | k — 1) are calculated as follows: 

1. The derivation of probabilities vector /xj“(.) = 
[^i(k) . . . n~(k)]' in classical IMM is replaced by 
the derivation of the predicted bba m)7(.) given by 

m fc (•) - M t ' m fe-r(-) (11) 

2. The derivation of probabilities fii\j(k — l|fc — 1) = 
P{Mi{k — l)\Mj(k),Z k ~ 1 } is replaced by the 
derivation of bba m k _ i|fc-i(-) thanks to the Gen- 
eralized Bayesian Theorem (GBT) [18]. More pre- 
cisely, 

m fc^i|fc_ 1 (.|M : ,(fc)) = 

[©mf 0(fc_1)xe(fe) (.|M i (fc - l))](.|Mj(fc)]l e ^ fc_1 -* 

(12) 

where -fr 0(/c — 1) x 0(/c) is the ballooning exten- 
sion [18] of the bba on the Cartesian product frame 
Q(k — 1) x 0(/c), and where j Q(k — 1) represents 
the marginalization operation of the bba on the 
frame 0(fc — 1). See [18], for details and examples. 

3. The derivation of the mixing probability Hi\j{k — 
1 | k — 1) = P{Mi(k — 1)| Mj(k), Z k ~ 1 } of classical 
IMM is replaced by the pignistic probability drawn 
from m fe _ 1 | fc _ 1 (.|Mj(/c)), that is: 

/i Aj {k - 1 | k - 1) = BetP{Mi{k - l)|M,-(Jfe), Z^ 1 } 

where BetP{.} is calculated with the transforma- 
tion (10) using m fe _ 1 |*,_ 1 (.| Mj(k)) given by (12). 

x 9 (fc — l\k — 1) and P 9 (k — l\k — 1) are calculated as in 
IMM Step 1. 

• Step 2: Same as IMM Step 2. 

• Step 3 (Mode bba update): The updated bba nifc(.) 

of modes is computed from the conjunctive combination 
of the predicted bba (.) given in (11) with observed 

bba’s 8 j = 1, 2, . . .r by 

m fe (.) = [m fc> i© . . . ©m fc>J .©m~_ 1 ](.) (13) 

where the observed bba’s for j = l,...,r are 

given 9 by [13]: 

m k ,j{M j (k)) =0 

m k,j(Mj(k)) = aj(l - RAj(k)) (14) 

mfc,j(0(fc)) = 1 — aj(l — RAj(k)) 

cXj is a discounting coefficient associated with the like- 
lihood of the mode Mj(k ) and R is a normalization 
constant. 

8 We mean that the bba is built from the likelihood 

A j(k) which depends on the mode Mj ( k ) and on the observation 
available z (k). 

9 This is Appriou’s model no. 1 in [1]. 



• Step 4 (Global estimation for output purpose): 
The global estimate x(/c|fc) and the covariance of 
estimation error P(/c|fc) are given as in step 4 of clas- 
sical IMM by taking Hj(k) = BetP{Mj(k) | Z fc } where 
BetP{Mj(k) \Z k } is the pignistic probability that the 
mode Mj is effective at time k. BetP{Mj(k)\Z k } 
is computed from the updated bba mfc(.) given by (13). 

A mistake in Step 1 of BIMM filter: The afore- 
mentioned Step 1 of BIMM algorithm described with an 
example in [13] is clearly incorrect because the deriva- 
tion of the predicted bba m]7(.) by (5) is wrong be- 
cause the sum of masses of focal elements is not equal 
to one. It is easy to verify from example in [13] 
when considering only two models, when taking pi = 
m.(Mi(k)\Mi{k — 1)) = 0.9, 1 — /5i = 0.1 = m^Mpk) U 
Mi(k)\Mi(k - 1)) and p 2 = m { M 2 (k)\M 2 (k - 1)) = 
0.89, 1 — p 2 = 0.11 = m^M 2 [k) \JM 2 {k)\M 2 {k — 1)) and 
taking the prior bba mfc_i(.) = [m(0) = 0 m(Mi(k — 
1)) = 0.45 m(M 2 (k — 1)) = 0.20 m(Mi(fc - 1) UM 2 (k - 
1)) =0.35]'. Applying the wrong formula (11), one gets 
precisely: 



'1 


0 


0 


0 ' 




0 ' 




0 




0 ' 


0 


0.9 


0 


0.1 




0.45 




0.4400 




0.44 


0 


0 


0.89 


0.11 




0.20 


— 


0.2165 


0.21 


0 


0 


0 


1 




0.35 




0.3500 




0.35 



M f mfc_i(.) mfc (.) Result in [13] 



One can see that the sum of components of m)7(.) 
equals 1.0065 !!! This mistake is not due to round- 
ing approximation of the result, but to a more serious 
mistake in the choice of the transition matrix M t . This 
mistake actually comes from the confusion in indices of 
the classical IMM transition matrix. It is easy to ver- 
ify that the correct transition matrix must be actually 
taken as the transpose of M t . Therefore, the correct 
derivation of m]7 (.) must be done by 

( 15 ) 

For the example 1 of [13], one will get correctly 



'10 0 O' 




0 ' 




0 


0 0.9 0 0 




0.45 




0.4050 


0 0 0.89 0 




0.20 




0.1780 


0 0.1 0.11 1 

V 




0.35 




0.4170 



— 1(0 m“(.) 



Remarks on BIMM filter: The BIMM is based 
on two 10 pillars: 1) the conjunctive rule of combina- 
tion, and 2) the pignistic transformation to approxi- 
mate a bba into a subjective probability measure be- 
cause. These two pillars are disputable because: 

10 Actually, Smets’ Generalized Bayesian Theorem (GBT) 
could be also considered as the third pillar of BIMM. 
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1. The efficiency of Smets’ rule for combining bba’s 
is very questionable in this belief-based extension 
of IMM because it has been already proved in [15], 
Vol. 3, and specially in sequential Target Type 
Tracking problem [7] that such rule doesn’t per- 
form well in general for mode change detection. 
Smets’ rule doesn’t respond to new information 
since very quickly all the mass of belief concen- 
trates on the empty set. See example in [15], Vol. 
3, Chap. 1, freely downloadable from the web and 
not reported here due to space limitation. 

2. The real interest and efficiency of the pignistic 
transformation is also disputable because there ex- 
ists other probabilistic transformations which per- 
form better than BetP in term of probabilistic in- 
formational content, in particular the DSrnP trans- 
formation developed in [15], Vol. 3, Chap 1 & 3 
and also in [8]. 

3. The justification for the use of Appriou’s model no. 
1 in step 3 of BIMM is missing and probably other 
(and maybe better) models could be developed to 
derive the updated bba mfc(.). This question has 
not been investigated in this paper and will be a 
source for future research. 

Interest of BIMM w.r.t. IMM: The potential ad- 
vantage of the belief-based IMM approach is to offer 
some robustness of the filter when replacing the strong 
constraint on the knowledge of probability of transitions 
7 Tij (usually based on ad-hoc assumptions on the mean 
sojourn time of the target in each mode) by a more flex- 
ible constraint on the transitions based on (very sim- 
ple and less specific) uncertain implication rules. With 
BIMM, one can also relax the knowledge of the prior 
probabilities of the modes by starting the tracking di- 
rectly with a vacuous belief prior of the modes. Of 
course, if one has good reasons to use a given prior 
of modes, this can be done easily in belief-based IMM 
approach which is also a nice features of such filter. 

4 PCR-BIMM algorithm 

To preserve the potential advantages of BIMM and 
to overcome its aforementionned problems, we propose 
to keep its general structure as a belief-based extension 
of classical IMM but we replace Smets’ rule by the more 
effective Proportional Conflict Redistribution rule no. 

5 (PCR5), or eventually the more simple PCR rule no. 

6 (PCR6), and to replace the pignistic transformation 
by the more effective DSmP transformation to estimate 
modes probabilities required in the IMM filter. We call 
this new algorithm, the PCR-BIMM filter. Before giv- 
ing the sketch of our PCR-BIMM filter, we just recall 
what are the PCR5 fusion rule and the DSmP trans- 
formation. All details, justifications with examples on 
PCR5 and DSmP can be found freely from the web in 
[15], Vols. 2 & 3 and will not be reported here. 



4.1 PCR5 and PCR6 fusion rules 

In DSmT (Dezert-Smarandache Theory) framework, 
the Proportional Conflict Redistribution Rule no. 5 
(PCR5) is used generally to combine bba’s. PCR5 
transfers the conflicting mass only to the elements in- 
volved in the conflict and proportionally to their indi- 
vidual masses, so that the specificity of the informa- 
tion is entirely preserved in this fusion process. Let 
TOi (.) and m 2 (.) be two independent 11 bba’s, then the 
PCR5 rule is defined as follows (see [15], Vol. 2 for 
full justification and examples): m.pcR 5 ($) = 0 and 
VX G 2 e \ {0} 



wpcffi(^) = mi(Xi)m 2 (X 2 )+ 

Xi,X 2 £2 e 

x 1 nx 2 =x 

r mi(A) 2 m 2 (A 2 ) m 2 (X) 2 m 1 (X 2 ) , , . 

Z - e [ rm(A) + m 2 (X 2 ) + m 2 (X) + m^X ,) 1 [ 

X 2 £2 
x 2 nx=0 



where all denominators in (16) are different from zero. 
If a denominator is zero, that fraction is discarded. Ad- 
ditional properties of PCR5 can be found in [9]. Exten- 
sion of PCR5 for combining qualitative bba’s can be 
found in [15], Vol. 2 & 3. All propositions/sets are 
in a canonical form. A variant of PCR5, called PCR6 
has been proposed by Martin and Osswald in [15], Vol. 
2, for combining s > 2 sources. The general formu- 
las for PCR5 and PCR6 rules are given in [15], Vol. 2 
also. PCR6 coincides with PCR5 when one combines 
two sources. The difference between PCR5 and PCR6 
lies in the way the proportional conflict redistribution 
is done as soon as three or more sources are involved 
in the fusion. For example, let’s consider three sources 
with bba’s mi(.), m 2 (.) and m 3(.), A fl B = 0 for the 
model of the frame 0, and mi (A) = 0.6, m 2 (B) = 0.3, 
m 3 {B) = 0.1. With PCR5 the partial conflicting mass 
m\{A)m 2 (B)m 3 {B) = 0.6 • 0.3 • 0.1 = 0.018 is redis- 
tributed back to A and B only with respect to the 
following proportions respectively: x^ CR5 = 0.01714 
and Xg CR5 = 0.00086 because the proportionalization 
requires 



~,PCRb 

X A 



TOi (A) 



y,PCR5 



m 2 (B)m 3 (B) 



m 1 (A)m 2 (B)m 3 (B) 
m.i(A) + m 2 (B)m 3 (B) 



that is 



thus 



n PCR5 



0.6 



„PCR5 



0.03 



0.018 



0.6 



x RCm = 0.60 • 0.02857 : 
x RCm = 0.03 • 0.02857 ; 



-0.03 
0.01714 
0.00086 



0.02857 



With the PCR6 fusion rule, the partial conflicting mass 
m\{A)m 2 {B)m 3 {B) — 0.6 • 0.3 • 0.1 = 0.018 is redis- 
tributed back to A and B only with respect to the fol- 
lowing proportions respectively: x RCR6 = 0.0108 and 

11 1. e. each source provides its bba independently of the other 
sources. 



325 




Advances and Applications of DSmT for Information Fusion. Collected Works. Volume 4 



„pcm 

■ n 

is done as follows: 



= 0.0072 because the PCR6 proportionalization 4.3 Sketch of P C R - B I M M 



„PCR6 
X A . 

'mi(A) 

that is 

„PCR6 

X A 

0.6 



r PCR6 
l B, 2 



y.PCRd 
l B, 3 



mi (A)m ,2 {B)m 3 ( B ) 



m 2 (S) m 3 {B) m.i(A) + m 2 (B) + m 3 (B) 



r PCR6 
l b, 2 

0.3 



~PCR6 

l B,3 

0.1 



0.018 



0.6 + 0.3 + 0.1 



= 0.018 



We briefly summarize the five steps of our PCR- 
BIMM filter. 

• Step 0 (Initialization at k = 0): Same as Step 0 of 
BIMM. 

• Step 1 (Interaction-mixing): Same as Step 1 of 
BIMM except that the predicted bba m)7(.) is com- 
puted by (15) instead of (11), that is 



thus 

! x pcr& =0.6-0.018 = 0.0108 

x B% m = °- 3 ' °- 018 = °- 0054 

x P. c 3 m = 0.1 • 0.018 = 0.0018 

and therefore with PCR6, one gets finally the following 
redistributions to A and B: 
x pcrg = 0 0108 

X PCR6 = X PCR6 + X PCR6 = q .0054 + 0.0018 = 0.0072 

From the implementation point of view, PCR6 is much 
more simple to implement than PCR5. For conve- 
nience, Matlab codes of PCR5 and PCR6 fusion rules 
can be found in [15, 16]. 

4.2 The DSmP transformation 

The DSmP probabilistic transformation is a seri- 
ous alternative to the classical pignistic transformation 
which allows to increase the probabilistic information 
content (PIC), i.e. to minimize the Shannon entropy, 
of the approximated subjective probability measure 
drawn from any bba. Justification and comparisons of 
DSmP(.) w.r.t. BetP(.) and to other transformations 
can be found in details in [8, 15], Vol. 3, Chap. 3. 
DSmP transformation is defined 12 by DSmP e (%) = 0 
and MX € 2 e \ {0} by 



m fc (.) = M(-m fc _i(.) (18) 

and the derivation of the mixing probability — 

l|fc — 1) = P{Mi(k — l)\Mj(k), Z fc-1 } of classical 
IMM is replaced by the DSmP probability drawn from 
m fc ^i| fc _i(.|M i7 -(fe)), that is: 

Hi\ j(k — 1] k — 1) = DSmP e (Mi(k — 1)| Mj(k), Z fc_1 ) 

where DSmP e (.) is calculated with the transformation 
(17) using m fc _ 1 | fc _ 1 (.| Mj(k)) given by (12). 

• Step 2: Same as IMM Step 2. 

• Step 3 (Mode bba update): The updated bba mfc(.) 
of modes is computed from the PCR.5 (or eventually 
PCR6) rule, denoted ®, of the predicted bba m)T_ 1 (.) 
given in (15) with bba’s mkj(-), j = 1, 2, . . .r by 

m fe (.) = [m M © . . . © nik,r © m“_ 1 ](.) (19) 

where the observed bba's m^ i j(.) for j = 1 ,...,r are 
given as in BIMM by (14). 

• Step 4 (Global estimation for output purpose): The 
global estimate it(k\k) and the covariance of estimation 
error P(fc|fc) are given as in step 4 of classical IMM 
by taking fij(k) = DSmP e {Mj{k) |Z fc } computed from 
the updated bba mjt(.) by (17). 



DSmP e (X) = 

re 2 e 



Y m ( z ) + e • c ( x n Y ) 



zcxcy 

C(Z) = 1 



E 



m(Y ) 

m(Z) + e • C(Y) 



ZCY 
C(Z) = 1 

(17) 

where C(X fl Y) and C(Y) denote the cardinals of the 
sets X fl Y and Y respectively; e > 0 is a small number 
which allows to reach a highest PIC value of the ap- 
proximation of m(.) into a subjective probability mea- 
sure. Usually e = 0, but in some particular degen- 
erate cases, when the DSmP e - o(.) values cannot be 
derived, the DSmP e> o values can however always be 
derived by choosing e as a very small positive number, 
say e = 1/1000 for example in order to be as close as 
we want to the highest value of the PIC. The smaller e, 
the better/bigger PIC value one gets. When e = 1 and 
when the masses of all elements Z having C(Z) = 1 are 
zero, DSmP e = i(.) = BetP(.). 



12 Here we work on classical power-set, but DSmP can be de- 
fined also for working with other fusion spaces, hyper-power sets 
or super-power sets if necessary. 



Remark: This preliminary version of PCR.-BIMM is 
perfectible because it still shares several points with 
BIMM 13 . In particular, the Step 3 of PCR-BIMM cal- 
culates, as in BIMM, m^j(.) with a model based on 
likelihoods A j(k) whose strong justification is missing. 
Further investigations will be done to improve this step 
3, as well as the Step 1 to get better performances of 
PCR-BIMM (if possible) in a future research. 

5 Simulation results 

In this section, we present the application of the 
PCR-BIMM to a ground target tracking problem. We 
consider a vehicule localized in (1000m, 5000m) in the 
cartesian referential (X, Y). We simulate a ground sen- 
sor located in (0,0) which is able to detect the moving 
target in range p and azimut 9. The gaussian mea- 
surement noise is supposed to be white and centered 
with the covariances a p = 20 m and erg = 0.008 rad. 
The sampling time is fixed to 2 seconds. For tracking 
the ground target we only consider two motion models. 

13 In particular, the GBT is still used in Step 2 of PCR-BIMM. 
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Estimated trajectory 




Figure 1: True target trajectory and estimated trajec- 
tories. 



detection and in the same time and get higher proba- 
bility when the target is moving in the same mode. 





A constant velocity motion model called CV 1, with a 
small noise acvi = Ito.s ” 2 and another constant ve- 
locity motion model called CV 2, with a bigger noise 
<jcv 2 = 4 m.s to palliate the target maneuver. The 
initial state for each IMM, BIMM 14 and PCR-BIMM is 
the true initial target state x(0). The transition Matrix 
P f is equal to : 



P f 



0.95 0.05 
0.05 0.95 



( 20 ) 



and the mass transition matrix M* for the BIMM and 
PCR-BIMM is same as in the paper [13]. The initial 
motion model mass is represented by the vacuous mass 
function. 

To compare the performances between the algorithms 
we used the root mean square error (RMSE) in loca- 
tion and velocity (figure 2) and the mean of the motion 
models probability obtained with 100 Monte-Carlo runs 
(figures 3, 4, 5). The first remark is, there is no signif- 
icant improvement by using the belief function in the 
IMM. In fact, the RMSE of the IMM, BIMM and PCR- 
BIMM are globally the same. However, we can observe 
a short difference of the PCR-BIMM error after the tar- 
get maneuvers between the time intervals [20, 30] and 
[40,50]. This observation carries along the second re- 
mark: the motion model transition duration is longer 
with the IMM (figure 3) and BIMM (figure 4) than the 
PCR-BIMM (figure 5). Then with the taken parame- 
ters for this simulation, the PCR-BIMM appears to be 
a good and fast detector of the motion models transi- 
tion. However, its computed motion models probability 
is inferior to the probability obtained with the IMM and 
BIMM. More investigations need to be done to see if it 
is possible (and how) to improve PCR-BIMM in order 
to preserve both the good performance of the maneuver 

14 Our BIMM implementations uses algorithm described in sec- 
tion 3 with (15) and additional normalization step mj,(.) in (7) 
since otherwise the BIMM algorithm doesn’t work at all due to 
the problem mentioned in section 2. 



Figure 2: Root Mean Square Error. 



IMM 




Figure 3: Motion Model Probability of the IMM. 



BIMM 




Figure 4: Motion Model Probability of the BIMM. 



PCR-BIMM 




Figure 5: Motion Model Probability of the PCR-BIMM 
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6 Conclusions 

In this paper, we have examined in details the re- 
cent BIMM algorithm and have corrected a mistake 
in it, and also identified some of its limitations. To 
palliate the problems of BIMM algorithm, we have de- 
veloped a more efficient belief-based algorithm, called 
PCR-BIMM, based on the Proportional Conflict Re- 
distribution fusion rule and on the DSmP probabilistic 
transformation to replace the conjunctive rule and the 
pignistic transformation used in BIMM. The derivation 
of the predicted bba of modes done incorrectly in BIMM 
is also fixed in our PCR-BIMM filter. The perfomances 
of PCR-BIMM with respect to the (corrected) BIMM 
and to the classical IMM have been evaluated from a 
simple maneuvering target tracking scenario through 
Monte-Carlo simulations. The results obtained in this 
paper show the ability of the PCR-BIMM to track ma- 
neuvering targets and also to improve the maneuver 
detection. It is important to note that such PCR- 
BIMM filter can be considered as more robust than 
IMM since PCR-BIMM requires less specific prior in- 
formation than IMM. Nevertheless, PCR-BIMM pro- 
vides globally the same RMS estimation errors perfor- 
mances as those obtained with the classical IMM which 
requires more specific prior information. Application of 
PCR-BIMM for tracking multiple maneuvering ground 
targets in a battlefield surveillance context is under in- 
vestigation and results will be published in forthcoming 
papers. 

References 

[1] A. Appriou, Classification based on multi-sensor 
uncertain data-fusion, Agardograplr 337 on Mul- 
tisensor Multitarget Data Fusion, Tracking and 
Identification Techniques for Guidance and Con- 
trol Applications, (in french), 1996. 

[2] Y. Bar-Shalom, K. C. Chang, H. A. P. Blom, 
Tracking a maneuvering target using input estima- 
tion vs. the interacting multiple model algorithm , 
IEEE Trans on AES, Vol. 25, pp. 296-300, March 
1989. 

[3] Y. Bar-Shalom, X.-R. Li, Estimation and Track- 
ing: Principles, Techniques, and Software, Artech 
House, Boston, MA, 1993. 

[4] Y. Bar-Shalom, X.-R. Li, Multitarget-Multisensor 
Tracking: Principles and Techniques, YBS Pub- 
lishing, 1995. 

[5] H. A. P. Blom, An efficient filter for abruptly 
changig systems , Proc. of the 23rd Conf. on Deci- 
sion and Control, Las vegas, NV, USA, Dec., 1984. 

[6] H. A. P. Blom, Y. Bar-Shalom, The Interacting 
multiple model algorithm for systems with marko- 



vian switching coefficients, IEEE Trans on AC, 
Vol. 33, No. 8, pp. 780-783, Aug. 1988. 

[7] J. Dezert, A. Tcliamova, F. Smarandache, P. Kon- 
stantinova, Target Type Tracking with PCR5 and 
Dempster’s Rides: A Comparative Analysis, Proc. 
of Fusion 2006 Int. Conf., Firenze, Italy, July 2006. 

[8] J. Dezert, F. Smarandache, A new probabilistic 
transformation of belief mass assignment, in Pro- 
ceedings of Fusion 2008 Conference, Cologne, Ger- 
many, July 2008. 

[9] J. Dezert, F. Smarandache, Non Bayesian condi- 
tioning and deconditioning, Int. Workshop on Be- 
lief Functions, Brest, France, April 2010. 

[10] X. R. Li, Y. Bar-Shalom, Multiple-model estima- 
tion with variable structure, IEEE Trans, on AC, 
Vol. 41, No. 4, pp. 478-493, April 1996. 

[11] X. R. Li, Engineer’s guide to variable-structure 
multiple-model estimation for tracking , in Y. Bar- 
Shalom and W. D. Blair, eds., Multitarget- 
Multisensor Tracking: Applications and Advances, 
Vol. 3, Chapter 10, Artech House, Boston, pp. 499- 
567, 2000. 

[12] X. R. Li, V. P. Jilkov, Survey of maneuvering 
target tracking. Part V: multiple-model methods, 
IEEE Trans, on AES, Vol. 41, No. 4, pp. 1255- 
1321, Oct. 2005. 

[13] G. Nassreddine, F. Abdallah, T. Denoeux, A state 
estimation method for multiple model systems us- 
ing belief function theory, Proc. of Information 
Fusion conf. (FUSION 09), Seattle, Washington, 
USA, July 2009. 

[14] G. Shafer, A mathematical theory of evidence, 
Princeton University Press, 1976. 

[15] F. Smarandache, J. Dezert, Advances and applica- 
tions of DSmT for information fusion (Collected 
works), Vols. 1-3, American Research Press, 2004- 
2009 . http://fs.gallup.unm.edu/DSmT.htm 

[16] F. Smarandache, J. Dezert, J.-M. Tacnet, Fusion 
of sources of evidence with different impoHances 
and reliabilities, (submitted to Fusion 2010 confer- 
ence) . 

[17] Ph. Smets, Constructing the pignistic probability 
function in a context of uncertainty, Uncertainty 
in AI, Vol. 5, pp. 29-39, 1990. 

[18] Ph. Smets, Belief Functions: The Disjunctive 

Rule of Combination and the Generalized Bayesian 
Theorem, Int. Journal of Approximate Reasoning, 
Vol. 9, pp. 1-35, August 1993. 



328 




